Quantum oscillations and high carrier mobility in the delafossite PdCo02 
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We present de Haas- van Alphen and resistivity data on single crystals of the delafossite PdCo02. 
At 295 K we measure an in-plane resistivity of 2.6 /ir2-cm, making PdCo02 the most conductive 
oxide known. The low-temperature in-plane resistivity has an activated rather than the usual T^ 
temperature dependence, suggesting a gapping of effective scattering that is consistent with phonon 
drag. Below 10 K, the transport mean free path is ^20 /xm, approximately 10^ lattice spacings and 
an astoundingly high value for flux-grown crystals. We discuss the origin of these properties in light 
of our data. 



PACS numbers: 

A number of discoveries, for example graphene and the 
superconductivity in the sodium cobaltate system, have 
focused attention on the properties of two-dimensional 
conduction on triangular lattices. One of the most promi- 
nent families of such materials is the delafossites ABO2, 
where the A atoms (Pt, Pd, Ag or Cu) are directly 
bonded to each other in triangular-lattice sheets, sepa- 
rated by BO2 layers. Possibilities for B include Cr, Co, 
Fe, Al and Ni pLl[2]. The delafossite family includes band 
insulators and semiconductors ^ ^S^, transparent semi- 
conductors [4], metallic magnets [5], candidate thermo- 
electrics [6^, and magnetic and magneto-electric insula- 
tors [7 . There are also nonmagnetic metals that present 
a combination of strong anisotropy, high carrier density 
and high mobility that is unique for bulk materials |3j [8j . 
In the long term, the delafossites may be a candidate fam- 
ily for multilayer technologies benefiting from their broad 
range of physical properties. High purity is a prerequi- 
site for such developments, and so the extremely high 
conductivity of nonmagnetic delafossite metals is worthy 
of investigation. Here, we concentrate on the most ex- 
tensively studied material, PdCo02. 

PdCo02 crystallizes in the space group R3m{D^^). 
The lattice parameters of the hexagonal unit cell at 298 K 
are a = = 2.830 A, and c = 17.743 A jOj. Single crystals 
have been grown by several groups pHTT] . The Co atoms 
have formal valence +3 and configuration 3d^ , and calcu- 
lations place the Fermi level between filled pseudo-t2^ 
and empty pseudo-e^ levels p!2UT4] . The Pd sheets there- 
fore dominate the conductivity, leading to high transport 
anisotropy. The Pd atoms have formal valence +1, in 
configuration 4:dP~^bs^. A key issue in understanding 
the physics of PdCo02 is the relative importance of the 
M and 5s states to the conduction. We will return to 



this in light of our data. 

PdCo02 is nonmagnetic, and so shows none of the frus- 
tration effects commonly associated with magnetism in 
triangular lattices. Band structure calculations [121 US] 
and angle-resolved photoemission spectroscopy (ARPES) 
data indicate that it has a broad, half- filled conduc- 
tion band, leading to a single, nearly two-dimensional 
Fermi surface (FS) with a rounded hexagonal cross sec- 
tion. This is different from graphene or graphite, since 
the Pd atoms form a simple triangular array rather than 
a honeycomb lattice. 

At room temperature, PdCo02 is the most conductive 
oxide known. To measure the in-plane resistivity pab, 
we cut three needle-like samples from a 6 /im-thick sin- 
gle crystal, of length ~ 400 fim and widths 39, 42 and 
120 fim. The measured resistances of all samples were 
consistent with pa6(295 K)= 2.6(2) jaQ-cm: two and four 
times, respectively, less resistive than other oxides known 
for their high conductivity, SrMoOs and ReOa [T6[ [17] . 
Among the elemental metals, only Cu, Ag and Au have 
lower resistivities (1.7, 1.6 and 2.2 /il]-cm, respectively). 
The high conductivity of PdCo02 is especially surpris- 
ing considering that the C0O2 layers do not contribute 
itinerant carriers: the carrier density is 30% that of Cu. 

To study the origin of this high conductivity, we com- 
bined a de Haas- van Alphen (dHvA) study, using a piezo- 
resistive torque method [18], with high-resolution resis- 
tivity measurements. The crystals were grown as de- 
scribed in Ref. [9]. They are clean: the resistivity samples 
have pa6(295 K)/pab{0) ~ 250. Torque measurements 
were carried out between 50 mK and 20 K, on a crys- 
tal ^ 120 pm across and 30 pm thick; example torque 
data are shown in Fig. [l] That the measurement current 
was not heating the samples was checked with simulta- 
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FIG. 1: Upper panel: torque curves against H, at the in- 
dicated angles between H and the c-axis, and magnetore- 
sistance of a resistivity sample. The ticks on this curve and 
the 4.6° torque curve label the cycle number of the difference 
frequency: the troughs in the resistance oscillations match 
up with points where the torque oscillation amplitude (and 
likely magnetic interaction) is largest. Lower panel: Fourier 
transforms over the field range 8.1-15 T. 



neous measurement of dHvA oscillations in Sr2Ru04, a 
renormalized Fermi liquid with well-known quasiparticle 
masses. Using a rotatable sample stage, the field H could 
be rotated by an angle 6 from the crystal c-axis towards 
the F symmetry point of the first Brillouin zone (as il- 
lustrated in Fig. [2|. The temperature dependence of the 
dHvA effect was studied at ^ = —6°, and the depen- 
dence at T = 0.7 K. 

The data in Fig. [l] are consistent with PdCo02 being a 
high purity, high carrier concentration, quasi-2D metal. 
At each angle we observe two large dHvA frequencies 
that differ by only ~3%, and that vary approximately 
as 1/ cos(^), confirming the overall cylindrical form of 
the FS. We also observe significant magnetic interaction 
(MI), an effect seen in the highest purity metals [19]. 
Most torque curves show a common characteristic of MI, 
a strong difference frequency amplitude, visible in Fig. [l] 
as an oscillation matching the beating of the two high 
frequencies. The difference frequency also appears in the 



magnetoresistance, ruling out torque interaction as its 
origin. 

Considerable insight into the physics of PdCo02 can 
be obtained from detailed analysis of the dHvA data. 
Following [20 and [21 , the corrugations on nearly- 
cylindrical Fermi surfaces can be expanded into harmonic 
components with amplitudes k^u'. 



/Cf(^, kz) = ^ k^r, cos /J.(j) I 

117/ ^ 



COS udkz 
sin udkz 



where kz is the z component of the wavevector, d = c/3 
the interlayer spacing, and the azimuthal angle, /cqo is 
the average radius of the FS, and the remaining k^jy are 
the amplitudes of harmonic components of the corruga- 
tion. The components considered in this paper are illus- 
trated in Fig. [2] The triangular lattice of the Pd sheets 
allows components with /i = 0, 6, 12, ... with an even {i.e. 
cosine) dependence on kz. In addition, because PdCo02 
lacks mirror symmetry along the c-axis, /i = 3,9,15,... 
are allowed with odd {i.e. sine) kz dependence. 

Since /C;^ -independent warpings cannot be obtained 
from dHvA frequencies, we set /ceo = 0.040(3) and 
ki2,o = 0.007(2) A-i to match the ARPES Fermi sur- 
face |10]. With these amplitudes fixed, we then allow 
^00 5 ^01 5 and k^i to vary to obtain the best match 
to the dHvA frequencies, setting all higher-order ampli- 
tudes to zero. The FS cross-sectional areas in the fit to 
the data were calculated by numerical integration j22] . 
The fit and the resulting FS are shown in Fig. |2] The 
data constrain the relative signs of /cqi, /co2 and ksi] koi 
is taken to be > to match electronic structure calcula- 
tions. 

Numerical integration of the volume enclosed by the 
FS yields a carrier density of 1.005(1) per formula unit, 
based on room-temperature values for a and b. If a and 
b shrink by ~0.25% at low temperatures then the carrier 
density becomes almost exactly one: the FS of PdCo02 
arises, to high precision, from a single half-filled band. 

As mentioned above, a key issue for understanding 
PdCo02 is the relative importance of the Pd M and 55 
orbitals in states near the Fermi level. The dHvA masses 
[m* = {h^ /27r)dA/dE^ where A is the FS cross-sectional 
area] provide the first piece of evidence that the 5s com- 
ponent is more important than has usually been assumed. 
Our data yield m* = 1.45(5)me and 1.53(3)me for the 
29.3 and 30.7 kT orbits, respectively (at = -6°) [23], 
corresponding to an electronic specific heat of approxi- 
mately 1 mJ/mol-K^, in fairly good agreement with a 
measured value of 1.28 mJ/mol-K^ |9l[23]. These masses 
would be very low for a (i-orbital-derived FS. For compar- 
ison, Pd metal (which has atomic spacing 2.75 A, just 3% 
less than in PdCo02) has a mainly 55-derived electron- 
like Fermi pocket, compensated by a 4(i-derived open hole 
surface, and the dHvA masses of these surfaces are ~ 2.0 
and 5-lOme, respectively [2l[25]. The PdCo02 dHvA 
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FIG. 2: Top: Observed dHvA frequencies F, obtained by Fourier transform over the field range 8.1 to 15 T (circles; darker 
shading indicates larger amplitude), and a fit (small points). Bottom: the separate harmonic components included in the fit, 
and their amplitudes, keo and ki2,o are obtained from ARPES data [10 , and /coo, ^oi, ^02 and /csi from the fit. koi gives the 
overall Yamaji-type form to F{0) and ko2 an asymmetry between the upper and lower frequencies, ksi gives an asymmetry 
between ^ > and ^ < 0, which is allowed by the lack of mirror symmetry in the crystal structure. The Fermi surface that 
results from summing these components is shown at right. 



masses are much closer to that of the Pd 55 than the Pd 
M surface. 

The form of the FS warping also emphasizes the role 
played by the 55 component. The dominant warping is 
/coi, suggesting interlayer transport by direct Pd-Pd over- 
lap, which would rely on the extended 55 component to 
cross the 6 A interlayer separation. Hopping via the Co 
atoms would allow transport through the more compact 
d orbitals, but connects interlayer next-nearest-neighbor 
Pd sites (through the path Pd-O-Co-O-Pd, illustrated in 
Fig. 5b of [15 ), leading to a large ksi. Published calcu- 
lations indeed show bumps on the FS with significant Co 
character that give a large k^i [ElIIS]- The dHvA data, 
however, show that koi is far larger than ksi. 

Hybridization between the Pd Ad^^ and 55 orbitals has 
been discussed in [3l |13j |26] . Reports on electronic struc- 
ture calculations have tended to mention the likelihood 
of a 55 component at Ef but emphasize the M contribu- 
tion p!2HT5] . However in a calculation of PdCo02 with- 
out the C0O2 spacer layers, i.e. bare Pd sheets, a 55 
manifold that straddles Ep is clearly visible (Fig. 9 of 
Ref. [Ej), and remains apparent in full PdCo02 band 
structures. 

To analyze our data and compare with previous cal- 
culation, we performed electronic structure calculations 
on PdCo02 using the general potential linearized aug- 
mented plane wave method as implemented in the 
WIEN2K package. We calculated 20,000 k points in the 
full Brillouin zone, and studied the effects of an on-site 
repulsion /7eff at the Co atoms. For /7eff = 0, we re- 



cover previously reported results, and find dHvA masses 
of around 1.9me, larger than the measured value. Intro- 
ducing a modest /7eff reduces the calculated masses to- 
wards the experimental values and simultaneously gives 
a ratio between ksi and koi that is in better agreement 
with experiment (Table 1). The high sensitivity of the 
in-plane masses, which arise from Pd-Pd bonding, to the 
Co site Coulomb repulsion appears surprising. A possible 
interpretation is that the Pd 55 component is extended 
enough to be sensitive to the charge distribtuion in the 
C0O2 layers; this should be the subject of further inves- 
tigation. 

Although conduction by 55 states might account for 
the scale of the conductivity of PdCo02 at room tem- 
perature, the low temperature resistivity has additional 
intriguing features. Measurement is challenging be- 
cause of the tiny resistance of typical samples. By us- 
ing transformers mounted at 1 K to provide low-noise 
(150 pV/a/Hz) passive amplification, we obtained the 
data shown in Fig. |3] The difference between pc and pab 



TABLE I: LDA+U results, and experimental data, [/eff is the 
effective on-site Coulomb repulsion at the Co sites, is the 
dHvA mass for the orbit in the kz — plane, and the 
kz = ±7r plane. 

mp rriz koi (A-i) ksi (A-i) ksi/koi 

Ues = eV l.QOme 1.82me 0.0039 0.0172 4.4 

Uefi = 5 eV 1.66 1.52 0.0181 0.0029 0.2 

experiment 1.53 1.45 0.0107(3) 0.0010(3) 0.1 
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FIG. 3: In-plane PdCo02 resistivity pab{T) (top), and c-axis 
resistivity pc(T) (bottom), with fits. At low T, pc has a de- 
pendence, while pah is nearly temperature-independent. The 
inset in the pah panel illustrates ku^ the minimum wavenum- 
ber for Umklapp scattering. 



is striking: while pc shows a dependence at low tem- 
peratures, pah becomes almost temperature-independent. 

Below 30 K, Pab{T) is much better described by an 
activated form than the usual electron-phonon form 
(Fig. [3|. This may be a consequence of phonon drag, in 
which the phonon population is pushed out of equilibrium 
by charge carrier flow. It is a well-known phenomenon in 
low-temperature thermopower, but has been seen in re- 
sistivity only in the alkali metals [27 . In clean materials 
momentum is transferred to the lattice through Umk- 
lapp processes. For the alkali metals and for in-plane 



conduction in PdCo02, the Fermi surface is closed, re- 
sulting in an activation energy ksTu = hcku, where ku 
is the minimum wavenumber for Umklapp processes (il- 
lustrated in Fig. [3| and c is the speed of sound. For 
T <^ Tu^ electron-phonon scattering transfers momen- 
tum between electrons and phonons but does not relax 
the momentum of the combined system, ku ~ 0.5 
in PdCo02. Fits to Pab{T) for two samples yield activa- 
tion temperatures of 165 (Fig.[3| and 168 K (not shown), 
giving c ~ 4400 m/s. This is in fair agreement with spe- 
cific heat measurements, where the magnitude of the 
term, 0.0619 mJ/mol-K^ [9], yields, for isotropic sound 
propagation, c = 3650 m/s. This analysis works equally 
well for potassium [28] . 

Fermi liquids usually have a low-temperature resis- 
tivity due to electron-electron scattering. Ac in Pc{T) = 
pc(0) + AcT'^ is 2 • 10-^ n-cm/K'^. Pab{T) shows a very 
small resistivity upturn at low T, and the upper limit 
we can place on Aab is 1.0 • 10~^^ Q-cm/K'^. This is 
not a particularly small value for metals with electronic 
specific heats comparable to PdCo02 [29], but it is sig- 
nificantly smaller than Ac scaled by the low-T resistive 
anisotropy {pc{0) / pab{0) = 1000). In the alkali metals 
K and Na the closed FS geometry is known to partially 
suppress electron-electron Umklapp processes, reducing 
A by about a factor of ten [30 . Such a suppression is 
also likely to apply to Aab of PdCo02. 

The combination of resistivity and dHvA data 
also shows that transport electron-phonon coupling in 
PdCo02 is very small over an extended temperature 
range. In Ref. [9j, Pab{T) of PdCo02 for < T < 500 K 
is fitted to a sum of the Bloch-Gruneisen form and an 
Einstein mode contribution. These fits return Debye 
and Einstein temperatures Qjj and 0^;, and transport 
electron-phonon couplings scaled by the the plasma fre- 
quency squared, Xd^/^p ^E^tr/^p [31]. The dHvA 
data here allow determination of the plasma frequency: 
^l,xx = ne^/meo = (7.2(1) -10^^ sec-^)^ Performing the 
same fit over the range 60 < T < 300 K, where the lower 
bound excludes the region strongly affected by phonon 
drag, we obtain = 329 K and O^; = 1106 K, in good 
agreement with [9]. A^^^tr and XE,tr come to 0.043 and 
0.023, respectively, which are exceptionally low values: 
Cu has XD,tr = 0.13 and Ca, a metal with a particu- 
larly low dHvA mass of ~ 0.6me, 0.05 [32 . These low 
couplings may also be due to the FS topology, through 
a partial suppression of Umklapp processes at tempera- 
tures T ~ T[/. Careful calculation would be required to 
test this hypothesis. 

Our measured residual resistivities pabiT 0), be- 
tween 0.007 and 0.008 pft-cm for two samples, are very 
low, corresponding to a transport mean free path of 
20 /im, or, in a standard interpretation, a defect density 
in the Pd layers of one in 10^ [33 . This is a remark- 
able value for flux-grown crystals. To set it in context, 
the fractional quantum Hall effect has recently been ob- 



5 



served in ZnO-based heterostructures, where, after years 
of research, the transport mean free path has been in- 
creased to ~1.5 /im [34 . This is considerably less than 
what we report here for a naturally two-dimensional ma- 
terial, motivating the study of delafossites as a natural 
material for multilayer oxide technologies. 

In summary, we have observed exceptionally high con- 
ductivity in PdCo02 and mapped its Fermi surface in 
detail using the de Haas- van Alphen effect. Our data 
suggest that the conductivity of PdCo02 is due primar- 
ily to Pd 5s overlap. For in-plane transport, the con- 
tributions to resistivity from electron-phonon, electron- 
electron and electron-impurity scattering are all anoma- 
lously low. The effect of the closed Fermi surface topol- 
ogy on Umklapp processes seems to play a role and is 
worthy of investigation. It is also possible, however, that 
scattering is suppressed by some non-standard mecha- 
nism not considered in the analysis we have presented. 
Our results provide motivation for further work on this 
fascinating material. 
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SUPPLEMENTAL MATERIAL 

Mass determination. The Lifshitz-Kosevich fits used to 
determine tiie cyclotron masses at 6 = —6° are siiown 
in Fig. [4] The difference frequency probably appears as 
a result of magnetic interaction. If data below 1.3 K 
are excluded from the 29.3 and 30.7 kT fits, out of 
concern that magnetic interaction might affect the 
observed amplitudes, then masses of 1.50 and 1.66me 
are obtained, respectively. 

Specific heat. The electronic specific heat is given by 
7 = Cei/T = 7rk%m* /3dh'^ x V, where V is the volume 
per mole. If m* is taken to be l.SOmg at all /c^, then 7 
comes to 1.03 mJ/mol-K^. 
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FIG. 5: Resistive anisotropy pc/ Pah against temperature. The 
error from sample geometrical factors is 10%. 



Resistive anisotropy. In the relaxation time approxima- 
tion the expected anisotropy of a corrugated cylindrical 
FS is given by: 



Pel Pah 



1 -1 



where is the Kronecker delta [20 and d the interlayer 
spacing. Using the /c^j, from Table 1 of the main 
text yields pc/ pah = 241(13). The observed resistive 
anisotropy is larger (Fig. [5|, especially below 25 K, 
where it rises to more than 1000. The suppression of 
Umklapp processes described in the main text may 
play a role in enhancing the anisotropy, and the larger 
anisotropy at low temperatures suggests that the defect 
scattering cross sections are themselves anisotropic. 
Overall, however, the scale of the resistive anisotropy of 
PdCo02 should be the subject of further investigation. 



Tight-binding model. In the text of the paper we 
describe how strong interlayer next-nearest-neighbor 
(nnn) coupling, as would result from Co-mediated 
hopping, would give a large ksi. For confirmation, we 
performed a tight-binding calculation. The transfer 
integrals included are in-plane nn and nnn coupling, 
tnn = 1 and tnnn^ ciud interlayer nn and nnn coupling, 
tz and tzz- tnnn = "0.23 and Ep = 0.22 approximately 
(although not uniquely) reproduce the half- filling and 
hexagonal shape of the Fermi surface. We calculate the 
FS for two sets of parameters, and obtain the harmonic 
components: 

(1) {tnn.tnnn.tz^tzz) = (1,-0.23,0.042,0.011): 

keo = 0.039, koi = 0.0106, ksi = 0.0011 A-^. 

(2) {tnn.tnnn.tz^tzz) = (1,-0.23,0,0.04): 

keo = 0.039, koi = 0.0082, ksi = 0.0103 A-^ 
In the second set, where tzz dominates tz^ ksi > /cqi, in 
disagreement with the experimental data. 
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FIG. 4: Fits of dHvA amplitudes at = —6°, calculated over 
the field range 12.86-15 T, to the Lifshitz-Kosevich form, for 
the indicated frequencies. The masses resulting from the fits 
are indicated in the figure. 



Dingle analysis. In attempting a Dingle analysis, we con- 
cluded that the dHvA amplitudes are probably limited by 
phase smearing rather than impurity scattering: to ob- 
serve a Dingle mean free path of 20 jam would require 
a field homogeneity better than one in 10^. Such high 
homogeneity would be difficult to achieve even in the ap- 
plied field, and PdCo02 can contain Co-based impurity 
phases that are likely to increase the internal field in- 
homogeneity [9[ [11]. Our dHvA amplitudes, after the 
effects of finite temperature and measurement time con- 
stant are removed, are shown in Fig. |6] Impurity scatter- 
ing should yield a straightforward exponential decay, and 
fits to the data yield Dingle mean free paths of 5500 and 
7000 A for the 29.3 and 30.7 kT frequencies, respectively. 
However the exponential fits do not fit very well: the 
29.3 kT amplitude shows pronounced minima at ^ 0.08 
and 0.11 T~^, and the 30.7 kT amplitude appears to 
have minima at similar locations. This double- minimum 
form persists over a wide range of 0, and so is not due 
to a small FS warping: it is the reason that many of 



the frequencies in Fig. [2] are double peaks with a ~60 T 
sphtting. Nor is it due to magnetic interaction: it per- 
sists to T ~ 3 K. Phase smearing is a possibihty, with the 
minima being resonances where destructive interference 
is maximum, however the origin of this possible phase 
smearing is unclear. 
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FIG. 6: Amplitudes of the 29.3 and 30.7 kT peaks against 
obtained by Fourier transforming over windows of width 
0.00167 T~^. Rt is the Lifshitz-Kosevich temperature atten- 
uation form, with T = 0.7 K, and Rr is the attenuation due 
to the measurement time constant. The solid lines are fits to 
the standard Dingle attenuation; F is the dHvA frequency in 
tesla, /mfp the Dingle mean free path in A, and H the field in 
tesla. 



